Calculation of the bootstrap current profile for the 
TJ-II stellarator 



J.L. Velasco^, K. Allmaier^, A. Lopez-Fraguas^, CD. 
Beidler'^, H. Maassberg'^, W. Kernbichler^, F. Castejon^, 
J. A. Jimenez^ 

1 Laboratorio Nacional de Fusion, Asociacion EURATOM-CIEMAT, Madrid, 
Spain 

^ Institut fiir Theorctische Physik - Computational Physics, Tcchnische 

Universitat Graz, Association EURATOM-OAW, Graz, Austria 

3 Max-Planck Institut fiir Plasmaphysik, IPP-EURATOM, Greifswald, Germany 

E-mail: joseluis . velascoOciemat . es 



Abstract. 

Calculations of the bootstrap current for the TJ-II stellarator are presented. 
DKES and NEO-MC codes are employed; the latter has allowed, for the first time, 
the precise computation of the bootstrap transport coefficient in the long mean 
free path regime of this device. The low error bars allow a precise convolution 
of the monoenergetic coefficients, which is confirmed by error analysis. The 
radial profile of the bootstrap current is presented for the first time for the 
100_44_64 configuration of TJ-II for three different coUisionality regimes. The 
bootstrap coefficient is then compared to that of other configurations of TJ-II 
regularly operated. The results show qualitative agreement with toroidal current 
measurements; precise comparison with real discharges is ongoing. 



1. Introduction 

The bootstrap current is a neoclassical effect triggered by the radial gradients of the 
density and temperature in the presence of an inhomogeneous magnetic field, see e.g. 
Refs. [1, 2] and references therein. 

The control of the bootstrap current, and thus of the total parallel current, 
may lead to the possibility of continuous operation in tokamak overdense plasmas. 
In stellarators, it can provide access to improved confinement regimes, by means of 
control of the rotational transform profile. On the other hand, the bootstrap current 
may perturb the desired magnetic configuration produced by the external coils. This 
is specially important for shearless devices, such as W7-X [3] or TJ-II [4]. 

Indeed, one of the main lines of research at the flexible heliac TJ-II is the relation 
between confinement and the magnetic configuration [5, 6, 7, 8, 9]. Ref. [10] offers a 
review of results at TJ-II and other stellarators. Therefore, an estimation of finite- 
j3 effects on the characteristics of the magnetic configuration such as the rotational 
transform, the magnetic shear, the position of rational surfaces or the magnetic well, 
is of great importance. These effects are not expected to be large, since TJ-II is a 
heliac-type stellarator [4] consequently showing a small Shafranov shift of the magnetic 
surfaces, but still this should be quantified. 
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Additionally, an Electron Bernstein Wave (EBW) heating system is being in- 
stalled at TJ-II [11]. Apart from its main purpose of heating overdense plasmas, it 
may be employed to drive current as well. Then, the total current (bootstrap plus 
the one driven by EBW) could be used to try to tailor the rotational transform pro- 
file [12, 13]. 

The calculation of the bootstrap current is a numerical challenge for non- 
axisymmetric devices, since the error estimates for computations in the long-mean- 
free-path {Imfp) regime of stellarators are usually large: DKES [14] has a convergence 
problem in the Imfp and the Monte Carlo methods have an unfavorable scaling of 
the computing time with the collisionality. Examples of similar calculations for non- 
axisymmetric devices may be found in Refs. [15, 16]. The dependence of the bootstrap 
current on the configuration has been explored in several devices such as W7-X [17] 
and Heliotron-J [18]. The former presents momentum correction techniques applied 
to the calculation of the bootstrap current profile for two different configurations of 
W7-X. In the latter, the monoenergetic bootstrap coefficient is studied as a function 
of the bumpiness of the configuration. Ref. [19] also contains a comprehensive study of 
the numerical results for the bootstrap coefficient, as well as the other monoenergetic 
transport coefficients, for several stellarators. 

Despite its interest, the bootstrap current profile at TJ-II is not known with 
precision: the problem of lack of accuracy in the Imfp regime is especially relevant 
for this device [20, 21, 22], which is characterized by a very complex magnetic 
configuration. Recently, some effort has been put in calculating with precision the 
bootstrap coefficient, see Ref. [19] and references therein. One of the codes in Ref. [19] 
is NEO-MC [23], which was explicitly developed in order to overcome the problem 
of poor convergence in the Imfp regime. It combines the standard 5f method (see 
Ref. [23] and references therein) with an algorithm employing constant particle weights 
and re-discretizations of the test particle distribution. Apart from the monoenergetic 
bootstrap coefficient, NEO-MC also calculates the radial diffusion coefficient, required 
for calculating the ambipolar electric field, and the parallel conductivity coefficient, 
required for the momentum correction techniques [17, 24]. 

In this work, we show that NEO-MC is able to provide, for the first time, 
calculations of the contribution of the Imfp regime to the bootstrap current of TJ- 
II with high accuracy. We present, for the 100_44_64 magnetic configuration of TJ-II, 
computations of NEO-MC of the monoenergetic bootstrap coefficient at several radial 
positions for a wide range of collisionalities and radial electric fields. This allows, 
by means of energy convolution and momentum correction technique, to estimate 
the profile of the bootstrap current. A Monte-Carlo technique for calculating error 
propagation allows us to state quantitatively that the results shown are of significance. 

TJ-II is a flexible Heliac whose magnetic configuration can be modified, both on 
a shot-to-shot basis [5, 6, 7, 8] or continuously [9]. The bootstrap current may change, 
even qualitatively, for different magnetic configurations of TJ-II [21]. We have thus 
explored part of the set of magnetic configurations of TJ-II with DKES calculations. 

This work is organized as follows. In Section 2, we describe the calculation of 
the bootstrap current: we present the equations, we outline the difficulties that may 
arise when solving them and we briefly describe the algorithm employed. Then, we 
particularize the discussion to the TJ-II stellarator: in Section 3, we show the re- 
sults. We first discuss the monoenergetic calculations, then the profile of the current 
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2. Basic theory and calculation method 

In order to estimate the bootstrap current, one has to solve the Drift Kinetic Equation 
(DKE) that describes the evolution of the particle distribution function [14]: 

d_ 

dt 

where Vg is the guiding center drift velocity, -0 is the toroidal magnetic flux through 
the local magnetic surface and Cc is the collision operator. Here / is the normalized 
perturbation of the distribution function / and it is defined through the local 
Maxwellian distribution function /m: 

/^/m-./^. (2) 

The neoclassical ordering [2] allows to neglect the time-dependence of /, as well 
as its radial variation. The equations become local, i.e., the radial position enters only 
as a parameter. If the E x B drift included in Vg is treated as incompressible [25], the 
neglect of energy diffusion in the Lorentz collision operator allows for monoenergetic 
solution of the equations: the kinetic energy can also be considered a parameter. The 
remaining equation is then three-dimensional. 

Ignoring momentum-correction for the moment, the fluxes may be expressed as 
linear combinations of the gradients of density and temperature and of the electric 
field. The following expressions, together with their derivation, may be found for 
instance in Ref. [19]. For each species b, the radial particle flux and the bootstrap 
current are: 

(rb)_ fldn E, 3 1 dTb\ , 1 dTb 

(Jh-B) . /Idn E, 3 1 drb\ , 1 dTb 



ZhenBo \n dr 2Th dr J ^^Tb dr 

As usual, Zb, Tb and n are the charge number, temperature and density of species b, 
e is the elemental charge and i?r is the radial electric field; b = e are electrons and 
b = i are protons. Bq is the (0,0) Fourier harmonic of the field in Boozer coordinates 
and the brackets denote flux-surface average. The thermal coefficients at each 
radial position are calculated by convolution with a Maxwellian distribution of the 
monoenergetic coefficients: 

q{r,n,Ti,Tc,E,)^^ dx' e--\'+'^'-'^+'--^ Di^ir,,,* , (5) 
V"' Jo 

where a; = v/v^^^ is the particle velocity normalized by the thermal velocity of species b ; 
Si^2 is the Kronecker delta. The monoenergetic coefficients depend on the radially 
local magnetic field strength, and are calculated and stored in a database for fixed 
values of the collisionality v* = vR/vl and the electric field parameter Q = E^/ (vBq). 
Here, u is the collision frequency (see Ref. [19] for its explicit form), R the major radius 
of the device and t the rotational transform. There is one independent simulation for 
each value of r, v* and O, from which we obtain the three independent monoener- 
getic transport coefficients Dn, D31 and (note that Dn = D12 = I?2i = -D22 and 
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D31 — D32 — Di:} — D23) ■ 

Following Eqs. (4) and (5), one needs to compute and convolute the monoenergetic 
bootstrap coefficient D31 in order to calculate the thermal coefficients L]^^ and L32 
which relate the parallel current to the thermodynamic forces. This convolution must 
be done for a given radial electric field i?r, that is obtained from the ambipolar 
condition (rc)(i?r) ~ (ri)(£'r)- Thus, according to Eq. (3), the radial diffusion 
coefficient Du must also be computed in order to estimate i^'i ^-nd L12, which relate 
the radial particle fluxes to the thermodynamic forces. 

Finally, the Lorentz collision operator does not conserve momentum. For this 
reason, momentum-correction techniques are included in the calculation, following 
method 2 in Ref. [17]. Then, the linear relation of Eqs. (3) and (4) is not correct 
anymore (but still allows the discussion of the qualitative behaviour of the bootstrap 
current in terms of the transport coefficients and the plasma gradients). For these 
techniques, one needs also to calculate the parallel conductivity coefficient D33. 

One must note that the local ansatz underlying this approach is only partially 
fulfilled at certain positions of TJ-II [26, 27]. This makes the diffusive picture of 
transport only approximately valid. Nevertheless, relaxation of the local ansatz would 
make the calculation of the bootstrap current profile impossible in terms of computing 
time. Also the monoenergetic picture [28] breaks down for very large electric fields. 

2.1. The Sf method 

Solution of Eq. (1) in order to calculate the bootstrap current is a numerical challenge. 
The reason is that one has to calculate an asymmetry in the particle distribution 
function caused by a combination of plasma radial gradients and particle trapping 
between local maxima of the magnetic field strength. Although the trapped particles 
do not provide a large current, the passing particle population equilibrate with them 
via collisions. Then, when it comes the time to calculate this effect by estimating 
Djj by means of a Montc-Carlo code, one has to be able to estimate the non- 
compensation of the current carried by co-passing particles and that carried by 
counter-passing particles. Note that the trapped-particle fraction is a radially local 
quantity determined by the variation of the magnetic field strength in the local fiux- 
surfacc labelled by r. Therefore, the dependence of the transport coefficients (and 
specifically of the bootstrap coefficient) on the trapped-particle fraction is implicitly 
accounted for in Dij{r, i^* of Eq.(5). 

Additionally, in stellarators, particle trapping may be toroidal or helical [1], and 
these two components may add one to the other or cancel depending on the magnetic 
configuration. 

Finally, in Sf Monte Carlo methods, the contribution of each test particle to 
the monoenergetic coefficient is proportional to its radial excursion from the original 
magnetic surface. Therefore, the contribution of the trapped particles consists of 
relatively large terms which numerically cancel out. This is the issue that is addressed 
with NEO-MC. 

NEO-MC [23] is an improved 6f Monte Carlo method developed for computation 
of mono-energetic neoclassical transport coefficients in stellarators. This method uses 
a re-discretization procedure and importance sampling in order to reduce the variance 
of these coefficients. As a result, the CPU time required for a given accuracy of the 
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computation of the bootstrap coefficient scales as the mean free path to the power of 
3/2. This method introduces no bias from filtering out particles. In addition, it allows 
simultaneous computation of the bootstrap coefficient and the diffusion coefficient. 
NEO-MC has been benchmarked with other Monte Carlo and field-line-tracing meth- 
ods and with DKES, for different magnetic configurations and for a range of values 
of the coUisionality and the electric field [19]. The results are in good agreement, and 
we make use of both NEO-MC and DKES in the following. 



2.2. Error estimate 

Although NEO-MC (and DKES for not too low coUisionalities) provide small error 
bars, the total current may be close to zero: since the helical and the toroidal 
contributions, as in every classical stellarator, have opposite signs, they may exactly 
cancel out at a given coUisionality. Therefore some error estimate is required in order 
to assess if the uncertainty in renders the calculated currents inaccurate. In this 
section, we show a method for doing so. 

Although Eq. (1) is linear, Vg includes the radial electric field, which we calculate 
by means of the ambipolar condition by root-finding (we select among multiple roots 
according to a thermodynamic condition, see e.g. [29]). Therefore L^{E-[) are employed 
to calculate E^- itself. In this non-linear calculation, error propagation from the 
monoenergetic coefficients to the bootstrap current profile may become cumbersome. 

We therefore undertake a Monte Carlo error propagation. We start from a 
database of monoenergetic transport coefficients with their corresponding error bars, 
which we have previously calculated with NEO-MC or DKES. We then take a sample 
of the database: at every radial position p, for every value of v* , and fi, we give 
numerical values to I^n, Z?3i and D33. We do so by generating, at every point in the 
(p,z^*,r2)-space, three independent Gaussian random numbers according to the average 
values and error bars of the coefficients stored at the database. Then, convolution, 
momentum-conserving and root-finding techniques provide the ambipolar electric field 
and the neoclassical fiuxes for this particular realization of the neoclassical database. 
Averaging over 50 of such samples is usually enough to present the bootstrap current 
profile with a meaningful variance. Here, p= y'l/'/V'o = r /a is the normalized radial 
coordinate, where V'o is the toroidal flux through the last closed magnetic surface. 

This procedure may also be straightforwardly extended to account for the uncer- 
tainty in plasma proflles, although, as we show below, we do not consider error bars 
in our density and temperatures in this work. 



3. Calculation and results 

3.1. Magnetic configuration 

TJ-II is a medium size = 4 heliac with bean-shaped plasma cross-section and 
strong helical variation of its magnetic axis. The magnetic configuration employed 
for these calculations is the so-called 100_44_64, the most often employed during the 
experimental campaign of TJ-II. The major radius of this configuration is i? = 1.504 m, 
the minor radius is a = 0.192 m and the volume-averaged magnetic field is 0.957 T. The 
iota profile is fairly flat, with t(0) = — 1.551 and ((a) = — 1.650. The vacuum equilibrium 




Figure 2. Variation of B/Bq along tlie magnetic field line passing through 8 = 
(/> = at p = 0.5 . 
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is used here. In Fig. 1 we plot the contours of constant magnetic field strength B at 
the flux-surface given by p = 0.5, where 6 and </> are the poloidal and toroidal angles in 
Boozer coordinates. The structure of deep localized minima and maxima is apparent. 
Fig. 2 show the variation of B along a field line starting at = = 0, which clearly 
cannot be described with a few Fourier modes. Hence the convergence problems of 
DKES. 

3.2. Monoenergetic coefficients 

We have calculated, for the configuration described in the previous subsection, the 
monoenergetic coefficient -D31, together with Du and D33 for a wide range of the 
collisionality and the electric field parameter. We calculate the dependence of the 
monoenergetic coefficient with i/* (between 3 x 10~^ and 3 x 10^) for fixed values 
of ri between and 1 x 10^^. We have repeated these calculations at several radial 
locations, from p = 0.08 to p = 0.98. 

Part of the results, at middle radial position {p = 0.5) are shown in Fig. 3. The 
monoenergetic coefficient is presented there normalized to that of an axisymmetric 
tokamak of the same aspect ratio in the banana regime: 

D-,,.D,,/Dt,. Dt,^,A,^J^^. (6) 

where e is the inverse aspect ratio and mb the particle mass. The results for D^i 
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are qualitatively similar to those observed for LHD, CHS, W7-AS and others [19]: for 
v* >1, the collision frequency v is larger than the bounce frequency, hence there are no 
trapped particles and no seeding mechanism for the bootstrap current. Note that the 
monoenergetic coefficient is already close to zero for lower collisionalities than = 
for the bounce frequency is smaller than the effective-detrapping collision frequency, 
which may be defined as v/e. For very low collisionalities, as the ratio between 
the electric field and the collisionality i/* increases, the coefficient converges to an 
asymptotic value depending on the configuration and the radial position [30, 31]. The 
predicted asymptotic limit has been calculated by means of field- line integration [32] . 

Low error bars are obtained in the whole range of collisionalities calculated while, 
for this configuration, DKES cannot provide accurate values for collisionalities lower 
than 3 x 10~^. Actually, these small error bars push the convergence with collisionless 
asymptote towards collisionalities even lower than 3 x 10~^. In the case of very 
small (especially zero) electric field parameter, this seems to happen at collisionalities 
that are completely out of reach for our computing resources. Nevertheless, for 
the collisionalities of the plasmas presented in Section 3.3, the data of Fig. 3 allow 
for a precise calculation of the bootstrap current. Only under reactor conditions 
(n ~ 3 X lO^Om-^, T ~ 20keV) would i^* - 5 x 10"^ — 10"'' correspond to thermal 
particles. 

Additionally, for some values of such as 3 x 10~^, the monoenergetic coefficients 
do not converge to the collisionless limit. This effect was already been observed in 
Ref. [19], and is still to be understood. Nevertheless, one has to keep in mind that 
for large values of the electric field, the starting hypothesis of the calculation are not 
fulfilled anymore [25, 27]: the kinetic energy is not an approximately conserved quan- 
tity. Fortunately, the electric field in this region of TJ-II is low enough, specially for 
NBI plasmas, and these data will not be required. 

Special attention must be paid to the interpolation and extrapolation of the 
monoenergetic database and to the calculation of the integral in Eq. (5). This 
integration is made by means of Gauss-Laguerre of order 64, and interpolation by 
3-point Lagrange, with v* and fl in logarithmic scale. In extrapolation, as in Ref. [33], 
we have made use of the asymptotic limit where available, but also tried extrapolations 
such as -D3;^(i/* < 3 X 10~^) = £'31(1'* = 3 x 10~^). We have made sure that none of 
these choices affects the final result. 

Once we have fully covered the relevant range of collisionalities with enough ac- 
curacy, the next step in the calculation of the bootstrap coefficient is to convolute 
the results shown in Section 3.2 with a Maxwellian distribution function. Then, the 
solution of a linear system of moment equations for which the coefficients are different 
energy moments of the mono-energetic transport coefficients [17] allows for the calcu- 
lation of the ambipolar electric field and then the bootstrap current. 



3. 3. Current profile 

The bootstrap current has been estimated to be of the order of 1 kA for low density 
plasmas of TJ-II [22, 34]. Since the magnetic field of TJ-II has a very broad Fourier 
spectrum, calculation in the Imfp regime is very complicated, and the results presented 
large error bars for low collisionalities [20, 21, 22]. It has been argued that, due to 
the exponential in Eq. (5), the unacceptably large error bars for v* ^ 10^^, well into 
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Figure 4. Plasma profiles for a low-density plasma of TJ-II. 

the Imfp might not prcckide the calculation of the bootstrap current [22], but this 
was not quantified as we propose in Section 2.2. These neoclassical calculations were 
compared to the measurements of the toroidal rotation velocity of some protons and 
impurities by means of passive emission spectroscopy [35]. Although there was some 
qualitative agreement, the results differed by a factor 2. 

Recently, lithium wall coating at TJ-II [36] has given access to low impurity 
concentration and high density at plasmas heated with Neutral Beam Injection. For 
these plasmas, in which the Imfp regime is less relevant, the bootstrap current should 
be easier to calculate. Nevertheless, being smaller in magnitude, the significance of 
the results must be checked. 

We calculate the bootstrap current profile for three regimes of TJ-II, correspond- 
ing to low, medium and high density. Figs. 4. 5 and 6 show plasma profiles char- 
acteristic of the three regimes. Low density plasmas at TJ-II have peaked electron 
temperature profiles, with central values of above 1 keV and hollow density profiles. 
As the density rises, both profiles become flatter and the electron temperature be- 
comes lower, of about 200 — 300 eV for high density plasmas heated by means of 
Neutral Beam Injection. The ion temperature is rather flat under all conditions, of 
the order of 100 eV. For these plasmas with lithium-coated walls, taking the effective 
charge Ze// = 1 is a reasonable approximation. Note that, since for thermal particles 
u* (X n/T'^ , higher density means higher collisionality, for the temperature decreases 
(for electrons) or barely changes (for ions). 
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Figure 7. Calculated bootstrap current profile for tfie low density plasma. The 
error bars are calculated following Section 2.2 



We show the resuhs in Figs. 7, 8 and 9, separating the contributions from electrons 
and ions. The latter may be compared to toroidal rotation measurements at TJ-II [35], 
available for low density plasmas. 

Some general features may be extracted from observation of the three figures. 
First of all, the current is mainly carried by the electrons. This is a consequence of 
the decoupling of the ion and electron temperatures, the latter being systematically 
higher, in TJ-II plasmas. Thus, the electrons are much less collisional and their 
contribution is much higher. Also, for the high-density cases (and the low density 
at outer positions), the plasma is in the ion root [37]: a negative radial electric field 
partially cancels out the ion channel while adding an additional contribution to the 
electron channel, see Eqs. (3) and (4). The exception are the inner radial positions 
in the low-density case, where the electron root is realized. There, the positive radial 
electric field drives a relatively high ion parallel flow (with low damping from the 
small trapped-particle friction) and the electron flow is coupled to it [17]. Also near 
the gradient zone of high density plasmas the temperatures are similar and the ion 
bootstrap current is not negligible. 

The maxima of the bootstrap current density follow those of the temperature 
gradient, since — 1.5 Ll^^ is systematically much higher than L^^. This reflects 
the overshoot in the bootstrap coefficient for collisionalities between 10~^ and 10~^, a 
usual feature in stellarators [19]. The higher the electron collisionality, the lower the 
amplitude of the maxima. 

There is a switch in the sign of the bootstrap current at medium radius, 
corresponding to a switch in the sign of the electron contribution. It may be 
traced-back as the change of sign in the bootstrap coefficient, corresponding to 
the cancellation of the stellarator-like and the tokamak-like contributions. This 
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Figure 9. Calculated bootstrap current profile for the high density plasma. 
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happens for a given coUisionality at each radial position (at medium radius, around 
V* ^= 3 X 10^^, sec Fig. 3). For higher densities, the switch in sign happens for larger 
minor radius. The consequence of this fact is that the total toroidal current is positive 
(tokamak-like) for low density plasmas and becomes negative (stellarator-like) for high 
densities. Integration of Figs. 7, 8 and 9 yield the following total currents, with the 
error bars in parentheses: 0.55(2) kA for the ECH plasma, 0.314(8) kA for the medium 
density NBI plasma and —0.739(3) kA for the high density plasma. This change of sign 
in the toroidal current has been observed in TJ-II operation, and the values obtained 
here are comparable to that measured by Rogowski coils [34]. 

The current in Fig. 7 presents a discontinuous behaviour around p ~ 0.69. It 
corresponds to the transition from electron to ion root: two stable roots exist from 
/9 = 0.63 to p = 0.71 and the electric field is evaluated by root-finding, so the transition 
layer is very narrow. If a differential equation for the electric field [29] were solved 
instead, the transition layer would be broader and Fig. 7 would become smoother. The 
qualitative results would be the same, with a total toroidal current ranging between 
0.5 and 0.75 kA. These results remain also of significance if error bars of 20% are 
allowed in the density and temperature profiles. 

A consequence of the high coUisionality of the ions is that their contribution to 
the bootstrap current suffers from very low error bars, while those of the electron con- 
tribution are much higher. Therefore, comparison of these results with measurements 
by Charge-Exchange recombination spectroscopy (CXRS) is to be done. Actually, 
much of the error bars in Figs. 7, 8, 9 comes from the uncertainty in D33 coefficient 
for low coUisionality, which is calculated with DKES and employed in the momentum- 
correction techniques. Although still the results are of significance, the calculation of 
D33 with NEO-MC is underway. 

In Fig. 10 we show an estimate of how the bootstrap current modifies the ro- 
tational transform profile. We also plot the rationals present in the plasma for the 
100_44_64 configuration. In the vacuum iota range, calculated with VMEC [38], sev- 
eral rationals (14/9, 11/7, 8/5 and 13/8) are present. Since the total toroidal current 
is small in absolute value for TJ-II, the most external rational barely moves. Inner 
rationals do move for ECH plasmas, which have plasma gradients for small minor ra- 
dius. Note that negative currents lead to reduced and hence unwind the rotational 
transform. 

The calculation presented is a linear approximation to a non-linear problem: this 
current is, as stated before, the result of the combination of the non-uniformity of the 
magnetic field in the flux surfaces and the density and temperature profiles. Since the 
calculated current is not zero, it will modify the vacuum magnetic field. The calcu- 
lation of the bootstrap current should be iterated in the new magnetic configuration, 
and the process should be repeated until one achieves convergence (i.e., the calculated 
current does not vary from one iteration to the next). This type of calculation has 
been made for LHD [39], which is a high-/3 device with relatively large bootstrap cur- 
rent and large Shafranov shift. It has also been performed for simplified models of 
NCSX and QHS [40] , in order to study the preservation of the quasi-symmetry. Since 
TJ-II lacks any of these quasi-symmetries, is designed to have small Shafranov shift, 
and the pressures and calculated currents are low, it is not necessary to iterate further. 
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Figure 10. Rotational transform profile corresponding to the vacuum 
configuration and the three plasma profiles. 



Configuration name 




V(m3) 


.(0) 


i{a) 


100_32_60 


1.087 


0.934 


-1.423 


-1.517 


100_38_62 


0.971 


1.031 


-1.492 


-1.593 


100_40_63 


0.960 


1.043 


-1.510 


-1.609 


100_42_63 


0.931 


1.079 


-1.534 


-1.630 


100_44_64 


0.962 


1.098 


-1.551 


-1.650 


100_46_65 


0.903 


1.092 


-1.575 


-1.676 


100_50_65 


0.962 


1.082 


-1.614 


-1.704 


100_55_67 


0.964 


1.073 


-1.659 


-1.739 



Table 1. Main parameters of the configuration scan 



3.4- Configuration dependence of the bootstrap current 

Finally, we make a scan in a relevant part the configuration space of TJ-II. A thorough 
study of the dependence of the parallel fluxes on the magnetic configuration is beyond 
the scope of this work; and also NEO-MC computations, although feasible, are quite 
time consuming. We thus try to describe the bootstrap current of the configuration 
in terms of its monoenergetic coefficient at p = 0.25 and p = 0.7 (where the electron 
temperature gradients are larger for low and high density plasmas respectively) for 
17 = 0. We calculate so for seven different configurations which are regularly explored 
in TJ-II operation. We show in Table. 1 the main global parameters: when switching 
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Figure 11. Main Fourier iiarmonics of tiic configurations 100_32_60, 100_44_64, 
100.55.67 



from one configuration to another, the volume is held approximately constant, while 
the rotational transform profile is raised without changing the shear. The neoclassical 
transport of these configurations is not expected to change much, since the main 
Fourier terms in the description of the magnetic field remain unchanged: Fig. 11 
shows the radial profile of the main terms in the Fourier decomposition in Boozer 
coordinates (see Ref. [41] for a more comprehensive discussion), given by 

oo oo 

B{p,(f>,0)/Boip)= Y.^-"-^rXpyos{me-Nn4,). (7) 

?l— — oo 771—0 

For the sake of clarity, only the standard 100_44_64 configuration and the extremes 
of the scan (100.32.60 and 100.55.67), arc presented in the plot. The two main 
Fourier components relevant for the bootstrap coefficient are the helical curvature 
bi^-i and the toroidal curvature 61,0- These components vary a few percents during 
the configuration scan. The next contribution may come from the toroidal mirror 
60,1 [15, 18]. This is the component of largest relative variation during the magnetic 
scan, even reversing its sign. Nevertheless, 60^1 is small compared with the two 
curvatures and therefore it does not have a large absolute impact on the bootstrap 
current, as we will see below. The relative variation in the remaining components of 
the Fourier spectrum is negligible. 

The calculation is made with DKES: up to 150 Legendre polynomials and up to 
2548 Fourier modes have been employed in the description of the distribution function 
and the largest 50 Fourier modes have been kept in the description of the magnetic 
field on every surface. We have made sure that increasing the number of modes does 




Figure 12. Bootstrap current coefficient at p = 0.25 as a function of the 
collisionality 



not modify the presented results within the error bars. 

In Figs. 12 and 13 we show the coUisionahty dependence of the bootstrap coeffi- 
cient for three different configurations at p = 0.25 and p = 0.7. For the sake of clarity, 
only configurations 100_32_60, 100_44_64 and 100_55_67 will be shown in the figures. 
Since the D^^ coefficient has already been discussed, we focus on the differences be- 
tween configurations. Although the calculations suffer from large error bars for low 
collisionalities, one may see that the normalized coefficients coincide for large colli- 
sionality and differ only for small collisionalities. The differences are relatively (but 
not absolutely) larger for small minor radius. At this region, the toroidal and helical 
curvatures are relatively smaller, and therefore differences in the toroidal mirror trap- 
ping show up in the bootstrap coefficient. 

A 1/t dependence of the bootstrap current is thus to be expected (since D31 oc 
D^i/l) for the medium and high-density plasmas of Figs. 5 and 6. 

For low-density plasmas, the situation is different for three reasons: first of all, the 
electron collisionality is significantly lower. Also the electron temperature gradient is 
large, see Fig. 4, so the contribution of the term proportional to L§2 Eq. (4) becomes 
larger. Finally, this gradient is large also for inner positions, where the differences in 
D^^ are relatively larger, see Fig. 12. The three facts tend to maximize the difference 
between configurations. 




Figure 13. Bootstrap current coefficient at p = 0.7 as a function of the 
collisionality. 



Finally, the sidebands not shown in Fig. 11 might yield different radial diffusion 
coefficients for different configurations. Although this situation is not likely to happen, 
the ambipolar equation would yield different radial electric fields for given plasma 
profiles, and thus different bootstrap current profiles. This calculation is left for future 
works. 

4. Conclusions 

Calculations of the bootstrap current for the TJ-II stellarator are presented. The Sf 
code NEO-MC, together with DKES, has allowed for the first time for the precise 
computation of the neoclassical bootstrap transport coefficient in the long mean free 
path regime of the vacuum configuration of TJ-II. The low error bars have allowed for 
a precise convolution of the monoenergetic coefficients. Since the calculation has been 
repeated at several radial positions, we are able for the very first time to calculate 
the radial profile of the bootstrap current of several plasmas of TJ-II with sufficient 
accuracy for comparison with experimental results. A Monte Carlo method for error 
estimation allows us to confirm the significance of the results. 

The current profile is well understood in terms of the plasma profiles. Since the 
electrons are within the Imfp and the plateau regime for all the plasmas considered, 
the gradient in the electron temperature provides the main contribution to the total 
current density. 
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The results show quahtative agreement with previous calculations and 
measurements at TJ-II; future work includes comparison with CXRS and Rogowski 
coils measurements for real TJ-II discharges. 

We have calculated the bootstrap coefficient for several configurations usually 
operated at TJ-II. The results suggest that the dependence of the bootstrap current 
on the magnetic configuration should be small, specially for NBI plasmas. A more 
comprehensive study is left for a future work. 

This calculation will allow to quantify more accurately finite-beta effects in the 
studies on the relation between confinement and magnetic configuration. Additionally, 
these results must be kept in mind when estimating the EBW current drive at TJ-II. 
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